Transient Nucleation near the Mean-Field Spinodal 
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Nucleation is considered near the pseudospinodal in a one-dimensional (j> 4 model with a non- 
conserved order parameter and long-range interactions. For a sufficiently large system or a system 
with slow relaxation to metastable equilibrium, there is a non-negligible probability of nucleation 
occurring before reaching metastable equilibrium. This process is referred to as transient nucleation. 
The critical droplet is defined to be the configuration of maximum likelihood that is dynamically 
balanced between the metastable and stable wells. Time-dependent droplet profiles and nucleation 
rates are derived, and theoretical results are compared to computer simulation. The analysis reveals 
a distribution of nucleation times with a distinct peak characteristic of a nonstationary nucleation 
rate. Under the quench conditions employed, transient critical droplets are more compact than the 
droplets found in metastable equilibrium simulations and theoretical predictions. 

PACS numbers: 05.70.Ln; 64.60.Qb; 82.20.Db 



I. INTRODUCTION 

Standard nucleation theory predicts a single 
metastable equilibrium nucleation rate, which gives 
rise to an exponential distribution of nucleation times. 
A key assumption in the application of the standard 
theory is that the system is in metastable equilibrium. 
In practice, any system that has undergone a quench 
relaxes into metastable equilibrium. During this relax- 
ation, the system has a finite probability of initiating 
a decay to the stable phase via an isolated droplet. 
This process, referred to as transient nucleation, is of 
both practical and theoretical interest. Experimen- 
tally, transient nucleation has been observed in the 
laser melting of thin silicon films [l[, crystallization 
in amorphous alloys 0, Q> an( i liquid crystals 
Theoretical work has concentrated on mean-field kinetic 
descriptions [1, 1, 0, B II E^JIlJllL Studies of transient 



have shown a 



nucleation in </> 4 models [bj 
distribution of nucleation times with a distinct peak and 
exponential tail. 

In a mean-field system, the spinodal defines the limit 
of metastability. The spinodal divides the phase-diagram 
into two regions, one region exhibits a metastable phase 
while the other does not. In a system with finite but 
long-range interactions, the analogous limit of metasta- 
bility is defined by the pseudospinodal. As the inter- 
action range increases, statistics at the pseudospinodal 
converge to the corresponding mean- field spinodal values. 
As the pseudospinodal is approached, the metastable sys- 
tem exhibits diverging susceptibility, correlation length, 
and correlation times that are characteristic of a criti- 
cal point [13, HI] ■ Nucleation near a pseudospinodal is 
structurally different from classical nucleation near the 
coexistence curve. In classical nucleation, droplets are 
compact fluctuations of the stable phase with a sharp 
interface separating the interior and exterior [ijj]. How- 
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ever, near the pseudospinodal, droplets do not exhibit 
a well-defined interface nor a stable-phase interior [20l |. 
Evidence for nucleation near the pseudospinodal is ob- 
served in deeply quenched liquids [111, [22| and in solid- 
solid phase transitions [23|, HH . Our objective is to study 
transient nucleation rates and transient critical-droplet 
profiles in a system with long-range interactions near the 
pseudospinodal. 

To make these ideas precise, we consider purely dis- 
sipative dynamics described by a Markovian Langevin 
equation that is a sum of a deterministic drift term 
and a stochastic noise term with zero mean. We de- 
fine the metastable well to be the set of configurations 
that would follow the (noiseless) deterministic drift to 
a stationary configuration that is not the global energy 
minimum. The metastable well boundary 23 consists of 
configurations that, upon random perturbation, would 
each drift to the metastable minimum with probabil- 
ity 1/2. A similar definition is employed by Roy et al. 
and others [2^, [26|, [27|- Transient critical droplets are 
characterized by the most probable configuration on the 
metastable well boundary at a given time t since the 
quench. 

We define the nucleation time as the latest time after 
the quench such that the system configuration was lo- 
cated on the metastable well boundary; the nucleating 
droplet is the corresponding system configuration. A sys- 
tem has nucleated when there is a negligible probability 
of returning to the metastable well. 

Because the nucleation rate is an extensive quan- 
tity, transience is best observed when the relaxation to 
metastable equilibrium is slow and the system size is 
large. Our primary results are that in magnetic field 
quenches at a temperature below the critical tempera- 
ture, the transient nucleation rate is adequately described 
by a quasi-equilibrium theory. Near the pseudospinodal, 
transient droplets are more compact than those that nu- 
cleate in metastable equilibrium. 

In Sec. II, we introduce the model and review the 
field theoretic treatment of nucleation near the mean- 
field spinodal, the basis for the subsequent analy- 
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sis. Section III considers the time-dependent likelihood 
of non-equilibrium states characterizing relaxation to 
metastable equilibrium. In Sec. IV, we find the time- 
dependent nucleation rate and transient critical droplet 
as a perturbation about the metastable equilibrium crit- 
ical droplet. In Sec. V, we compare our theoretical treat- 
ment to results from computer simulation. Section VI 
interprets and compares our results to previous work. 



II. METASTABILITY NEAR THE MEAN-FIELD 
SPINODAL 

We consider a one-dimensional ferromagnetic system 
with long-range interactions prepared in equilibrium with 
an initially negative external field, hinmai < 0. At time 
t = 0, the system is quenched and the external magnetic 
field h is set to a positive value. The system then evolves 
toward metastable equilibrium. At some point after the 
quench, it will decay to the stable phase via nucleation. 

The potential of the system is given by 



V{4>) = e<f> 2 + - hep, 



(1) 



where e < and u > 0. For systems of interaction range 
R, the Ginzburg-Landau- Wilson Hamiltonian is 



effect on nucleation. Therefore we drop the quartic term 
in the subsequent analysis. Within this cubic approxi- 
mation, the location of the potential minimum is 




< 0. 



(5) 



We introduce the shifted field, ip(x) = 4>(x)—6 m i n , so that 
the minimum of the mean-field potential occurs precisely 
at ip = 0. We rewrite the truncated potential in terms of 
the shifted field and drop an overall constant to obtain 



vw = - 



^6h sp A h ^ 2 



20; 



sp 



I^ 3 = ^ 2 -^ 3 , (6) 



where the variables a = — ^/6/i sp A/j/ 2d sp > and 
b = a 2 1 3A/j > are introduced to simplify the notation. 
Near the mean-field spinodal, the evolution of the field 
in the neighborhood of the metastable well is given by a 
Langevin equation, 



dip 
~dt 



R 



d 2 4> 
dx 2 



2aip + Zbijj 2 + rj. 



(7) 



We define the stationary configurations by setting the 
time derivative and the noise term to zero, 



H\6\= / dx 



R 2 
2 



dx 



(2) 



We assume purely dissipative dynamics, given by the 
time-dependent Ginzburg-Landau equation, 



5H 



dt 



7/ = R 



dx 2 



2e6 - Au6 3 + h + i ll (3) 



where rj = r)(x, t) is a zero-mean white noise with 
(r](x,t)r](x',t')) = 2/3- 1 5(x - x')S{t - t'), and (3 is 
the inverse temperature. The mean-field metastable 
well vanishes when the applied magnetic field is the 
spinodal field, defined by h 2 p — 8|e| 3 /27u. The cor- 
responding mean-field spinodal magnetization satisfies 
6 2 p = \e\/ 6u. These quantities are obtained from set- 
ting dV/dS = d 2 V/dd 2 = 0. 

Under the given quench conditions, the spinodal field is 
positive, h sp > h > 0, and the mean-field spinodal mag- 

< 0. For convenience, we in- 



a measure of the depth of the 



nctization is negative, 
troduce = h sp — h, 
metastable well. 

We expand the potential about the spinodal magneti- 
zation, 6 spi retaining terms up to third order [20L [2a|. 



V = 4> sv \ A h 



l K 



sp 



''sp 



(4) 



2 4 2 P 



■ o((0 - <M 4 ). 



We concentrate on nucleation in the neighborhood of the 
metastable well. Near the spinodal, the details of the po- 
tential for large positive magnetizations has a negligible 



= i? 2 ^!- 2a ^ + 3W> 2 . 
dx 2 



(8) 



Equation (|8|) has two spatially uniform solutions, the sta- 
ble configuration, 4>min(x) = 0, and an unstable configu- 
ration. There exists another non-uniform stationary con- 
figuration ips{x) that corresponds to the pro file of the 
critical droplet in metastable equilibrium |l9l . |20| . Be- 
cause the droplet can appear anywhere in the system, 
we fix the center of the droplet at the origin. A droplet 
must be symmetric, ips{x) = ^si^x), and cannot ex- 
hibit a sharp peak at its center, ip' s (0) = 0. At large 
distances, the droplet profile approaches the metastable 
background, \\m x ^± 00 ips(x) — 0. Solving Eq. ([8]), the 
non-uniform solution is 



^As = 7 cosh 
o 



a 
2~R2 



(9) 



This stationary configuration is unstable: assuming 
noiseless dynamics, a random perturbation about V's 
would cause the system to return to its metastable well 
with probability 1/2. Linear stability analysis about the 
non-uniform solution reveals a single negative eigenvalue, 
— 5a/2. The corresponding eigenvector gives the initial 
direction of growth (or decay) of the metastable equilib- 
rium critical droplet [19j, |28|, |29( , 



ipj\f = cosh x 



a 
2~R2 



(10) 



The amplitude of the growth eigenvector is concentrated 
at its center; a spinodal critical droplet grows by first 
filling its center [29| . 
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III. CHARACTERIZING RELAXATION TO 
METASTABLE EQUILIBRIUM 

In this section, we estimate the time-dependent like- 
lihood of field configurations, which is required to find 
transient critical-droplet profiles [3(1 HH, H2] • After the 
quench, the time evolution of the probability is given by 
the functional Fokker-Planck equation. We sidestep the 
formal solution and use a quasi-static mean-field theory. 
We assume at a time t since the quench, the likelihood 
of a configuration in the metastable well can be approx- 
imated by 



P[ip\t] k Cexpl -f3 / dx 



(11) 



where C is an undetermined constant, and F t = F t {jp) 
is a time-dependent polynomial in the field. At t = 0, 
Eq. (TTJ) should agree with the equilibrium distribution 
prior to the quench. At long times, we must recover the 
metastable potential, 



lim F t = V(i/j) = atp 2 - b^ 3 



(12) 



Because the quench only changes the external magnetic 
field, we assume that F t differs from the metastable equi- 
librium potential given in Eq. ([6]) by a time-dependent 
effective external field, J(t), 



F t = ai\} 2 - fo/> 3 + [h - J(t)]tp. 



(13) 

In the long time limit, Eq. (fT2|) implies that J(t) — ► h. 
We expect that the quartic term can still be neglected 
since this approximation will be used when t is close to 
t eq , when the system is in metastable equilibrium. 

We choose J(t) so that the density given in Eq. (fTTj) 
is consistent with the time evolution of the mean-field 
magnetization. In the mean-field limit, all fluctuations 
are suppressed and the system is described by a single 
scalar order parameter, ip(x,t) — » mt. Within the spin- 
odal approximation, the time evolution of the mean-field 
system is governed by the ordinary differential equation, 



d 
~dt 



—mt = —2am t + 3bm 2 . 



(14) 



The time evolution of the mean-field magnetization is 

2amo 



3&mo + (2a — 36mo) exp(2at) 



(15) 



where the initial condition mo denotes the prequench 
shifted mean-field magnetization. In order that Eq. (TTTj) 
reproduce the mean-field dynamics, Ft evaluated at m t 
must be a minimum, 



dF)_ 

dip 



= 0. 



(16) 



ip—mt 



The expression for the time-dependent effective field is 




FIG. 1: The mean transient magnetization prior to nucleation 
(solid line) obtained from the solution of Eq. (|3} compared to 
the average obtained from simulating Eq. (|19p in metastable 
equilibrium (squares). Both curves rise above the metastable 
mean-field minimum (broken line) given by Eq. (JT5J ■ 



J(t) = 2arrit — Sbm 2 + h. 



(17) 



Formally, Eqs. (jlip . (j 13[) . and (|17p define the first term 
of a large- R asymptotic expansion of the solution to the 
functional Fokker-Planck equation in the metastable well 
near the mean-field spinodal. 

Both the location of the minimum of F t and its cur- 
vature change as a function of t. For convenience, we 
introduce Aj = h sp — J(t) and 

A t = ~F t "(mt) = ay/Aj/A h , (18) 

a measure of the curvature of F t about its minimum. 
At long times, the mean-field magnetization approaches 
its metastable equilibrium value, mt — * 0, which implies 
that At — > a in the same limit. 

For a fixed time t, we define a corresponding fictitious 
Langevin equation, 

?*£A = R ^-F lm + *,,.,, (19) 

where the coordinate s is distinct from the time t since 
the quench. Once in metastable equilibrium, Eq. (|19p 
samples states from Eq. (jlip . We note that the states 
sampled by Eq. (fH)]) must be restricted to those within 
the metastable well of the complete dynamics of Eq. ([3]). 
For example, configurations that would nucleate Eq. (fH))) 
must be discarded because they are not in the metastable 
well of the underlying dynamics. With this restriction, 
we use Eq. (fl9|) to compare the statistics of Eq. (fTTj) 
at a fixed time t to those generated from simulation of 
the complete dynamics in Eq. (J3]) . Figure Q] compares 
the ensemble-averaged transient magnetization obtained 
from Eq. ([3]) to the metastable equilibrium magnetization 
sampled using Eq. (|19p for various t. 

To analytically estimate the transient magnetization 
of Eq. (fTTj) , we compute the stationary spatial average of 
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Qi 

FIG. 2: Schematic equipotential curves for a model system 
with two arbitrary dynamical variables, Qi and Qi- The 
broken line illustrates the linearization of the metastable 
well boundary. 



0.40 
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FIG. 3: The theoretical distribution of nucleation times 
computed from Eq. (|35[) (solid line) compared to the nor- 
malized histogram of nucleation times obtained from di- 
rect simulation of Eq. (J3J). 



Eq. CH), 



= -2a(ip) + 3b(ip 2 ) + J(t) - h, 



(20) 



which relates the first and second moments. The second 
moment can be approximated by the second moment of 
the Gaussian theory obtained by setting 6 = 0, 



1 



2(3R (2A t )!/2 



(21) 



Equations (|20j) and (|21j) yield an estimate of the magne- 
tization as a function of A t , defined in Eq. |TJ 



(ip) w m t + G t 1 > 



where, 



G t = ^RAha- 1 y/2A t > G c 



(22) 



(23) 



This estimate of the magnetization also agrees with the 
average magnetization obtained from simulation, shown 
in Fig. [TJ 



IV. TRANSIENT NUCLEATION 

At a time t after the quench, the transient critical 
droplet ipc(x,t) is determined by maximizing the proba- 
bility functional in Eq. (fTTj) constrained to the metastable 
well boundary 25, 



ipc(x,t) — maxP[V>|i]. 



(24) 



Note that the dynamically defined stationary config- 
uration ips given in Eq. lies on the boundary. Fur- 
thermore, we assume that our transient critical droplets 



resemble the dynamical fixed point, ips- In this limit, we 
approximate the metastable well boundary by a hyper- 
plane (25^ that is normal to the growth eigenvector ipw, 
given in Eq. (fTU)) . Figure [5] gives a schematic of the ap- 
proximation for a two-dimensional system. Within this 
approximation, a boundary configuration ip<g satisfies 



dx (ip<& - i/Jstytf = 0. 



(25) 



To find the configuration of maximum likelihood that 
satisfies the constraint, we introduce the Lagrange mul- 
tiplier At and extremize the functional, 



dx 



I? 
2 



dip 
dx 



+ F t (ip) - X t ip ipx 



(26) 



Centered at zero, the time-dependent critical-droplet 
profile solves the Euler-Lagrange equation, 



R 



d 2 ^ c 
dx 2 



= 2aip c - Zbipc - J{t) +h- \ t ipM, (27) 



where A t is chosen so that ipc satisfies the constraint in 
Eq. ([23)1 . For large t, the transient solution reduces to 
the dynamical fixed point, \im t — ,00 ">Pc = V's- At earlier 
times, the equation and constraint can be solved numer- 
ically. Because we have assumed the transient critical 
droplet resembles the fixed point, we expect At <C 1- 
This suggests a solution as a perturbation in At- First, 
we take At = and find the analog to Eq. |9]), 



ipc(x,t) = m t 




(28) 



Generally, the perturbation series is 

ip c (x,t) = ip° e (x,t) - At y(x,t) + 0(X 2 t ). (29) 
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FIG. 4: Ensemble-averaged nucleating-droplet magneti- 
zation profiles from simulation at t = 4.5 (solid line) and 
t = 7.5 (broken line). 



FIG. 5: Theoretical critical-droplet magnetization pro- 
files from Eq. ((29} computed at t — 4.5 (solid line) and 
t = 7.5 (broken line). 



The resulting differential equation for y(x) is 

R2 d^ = 2 v( x ) A t [1 " 3cosrT 2 (xVA/2i? 2 )] 
-cosh~ 3 (xVa/2i? 2 ). 



(30) 



We retain terms to first order in At, and substitute 
Eq. (|2"9")l into the constraint to find the following equation 
for the multiplier At, 



A, 



/ dx t) - ips) ijw 
J dx y(x) i/jm 



(31) 



The effect of the constraint is to decrease the amplitude 
of the transient critical droplet. 

With the transient droplet profiles characterized, we 
estimate the time-dependent nucleation rate. Given a 
time t before any nucleation event, there is unit proba- 
bility of being in the metastable well. With this normal- 
ization, the nucleation rate is proportional to the proba- 
bility to realize the transient critical droplet. We approx- 
imate this normalized probability by the relative likeli- 
hood of the transient critical droplet over the most likely 
metastable configuration, the uniform field with magne- 
tization rat-, 



r(t) 



P[Mx,t) 1 1] 

P[tp — m t | t ] 



(32) 



This rate estimate is simplified by defining the shifted 
effective potential, 



AF t (V0 = Ft(1>) ~ FtN 

= At^-mtf -36(^-m t ) 3 . 

With this definition, the transient nucleation rate is 



(33) 



Tf = T cxp <-/3 dx 



i? 2 /<M>_ 
~2~ Vdx 



AF t 



(34) 



where T is a prefactor that depends on the details of 
the system 33]. For the quench considered here, T(t) is 
a strictly increasing function of time. The experimentally 
accessible nucleation time distribution is given by 



p(t) = T{t) exp 



dt T(t) 



(35) 



V. NUMERICAL CALCULATIONS 

We simulated Eq. ([3]) for a periodic system of length 
L = 10 4 where e = -5/9, u = 1/4, R = 10, and /3 = 10. 
Prior to the quench, the system is prepared in equilibrium 
with an external field hinmai = —1.40 and shifted magne- 
tization mo = —0.727. After the quench, the applied field 
is set to h — 0.430. The distance to the mean-field spin- 
odal is Ah = 0.021, and the location of the metastable 
minimum is given by </> m in = —0.715. The corresponding 
parameters in the cubic approximation are a — 0.195 and 
b = 0.609. 

Throughout the run, the configuration of the system 
is saved periodically. Once the system has nucleated, we 
search for the nucleating droplet and nucleation time. We 
load a saved configuration and find the latest configura- 
tion that, upon perturbation, drifts to the stable phase 
with probability 1/2 [26j. We do not consider runs where 
multiple droplets appear at the nucleation time. The nu- 
cleation times were binned and compared to the theoret- 
ical results (see Fig. [3]) where the free parameter To was 
chosen to produce the best fit. Ensemble-averaged nucle- 
ating droplet profiles for two different times are plotted 
in Fig. 0] The corresponding theoretical critical droplets 
obtained from Eq. (|29|) are plotted in Fig. O 

Figures [3] and \5\ differ slightly because they repre- 
sent different statistical quantities. The cubic term in 
Eq. pip give rise to a skew in the distribution. Conse- 
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quently, average and extremal quantities do not agree. 
The same effect is demonstrated in Fig. Q]and Eq. (f2"2")l . 
where the transient system magnetization is found to dif- 
fer from the mean-field magnetization. For distributions 
that are sufficiently peaked, we expect the discrepancies 
to vanish [l^, l20l |. This occurs when the metastable- 
equilibrium Ginzburg parameter Goo given in Eq. (|23p is 
large. In the case considered here, the Ginzburg param- 
eter is finite, Goo = 26, resulting in the contrast between 
the averaged configurations and the extremal configura- 
tions. However, a large G M suppresses the nucleation 
rate, T ~ Lexp(— Goo), and obscures transience in the 
nucleating droplets. To preserve the transient-nucleation 
regime while increasing the Ginzburg parameter, the sys- 
tem size L must scale exponentially with Goo • 

VI. SUMMARY AND DISCUSSION 

We considered transient nucleation in a long-range one- 
dimensional <f> model with dissipative dynamics. We de- 
fined the metastable well boundary as the set of con- 
figurations balanced between the metastable and stable 
phases. We estimated the time-dependent likelihood of 
system configurations. We defined the transient criti- 
cal droplet as the most likely configuration constrained 
to lie on the metastable well boundary. We computed 
the likelihood of the transient critical droplet at vari- 
ous times. Our results explain the nonstationarity of nu- 
cleation rates reported in studies of the time-dependent 
Ginzburg Landau equation. 



As presented here, the theory is applicable to systems 
near metastable equilibrium where the system magneti- 
zation evolves slowly compared to the relaxation of other 
dynamical variables, as is the case near a spinodal. 

The theory provides a qualitative picture of the mea- 
sured ensemble-averaged profiles without any free param- 
eters. The analysis produces a distribution of nucleation 
times consistent with the simulation results with a single 
free parameter. 

Our analysis reduces to the earlier work [2(| at long 
times, when the system has relaxed to metastable equi- 
librium. We found that transient droplets decay to 
the background magnetization in the system at the 
time of nucleation. Before metastable equilibrium, the 
background magnetization acts as an anchor: transient 
droplets must have greater amplitude, which in turn sup- 
presses their rate of formation. Furthermore, this sug- 
gests that in comparable experiments, transient effects 
result in configurations that are more compact than pre- 
dicted by the metastable equilibrium analysis. Moreover, 
in systems with a stable crystalline phase, these results 
imply transience may determine the symmetry of critical 
droplets when nucleating near a spinodal. 
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